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Q_i| Abstract 

We show how the generalised (or skewed) parton distributions of the proton, H (x, £; kf,^ 2 ) 
unintegrated over the partonic transverse momenta, can be calculated from the known 
conventional parton distributions, q(x, fj?) and g(x, /j, 2 ), for small values of the skewedness 
parameter £. We demonstrate the procedure by numerically evaluating the skewed un- 
integrated gluon. We also provide a simple approximate phenomenological form of the 
distribution, which may be used to make more rapid predictions of observables. 



1 Introduction 

Generalised (or skewed) unintegrated parton distributions are essential ingredients in the cal- 
culation of many high energy processes initiated by protons. Examples of such processes are 
virtual photon Compton scattering (7*^ — >■ jp) ; diffractive heavy photon dissociation processes, 
like diffractive vector particle electroproduction (e.g. j*p — > pp, J/ipp, Tp, Zp) or dijet produc- 
tion (7*^ — > jj + p) or heavy quark production (j*p — > QQ + p); double diffractive central 
particle production (e.g. Higgs or dijet production, pp — » p + H + p, p + dijet + p). In these 
examples, the + sign is used to indicate the presence of a rapidity gap. 
We use the 'symmetrized' distributions introduced by Ji [ffll-0 



H a {x,q; k t , k t ; /i z 



(1) 



where a = q,g, with support — 1 < x < 1. They depend on the momentum fractions x\$ = 
x ± g^ carried by the emitted and absorbed partons, see Fig. 1. Since, in general, X\ 7^ x 2 , 
we speak of skewed or generalised distributions. We use the term unintegrated to indicate that 
the transverse momenta k t , k' t of the partons have not been integrated over. The bar on H 
is simply to denote that ([!]) is the unintegrated distribution. We reserve H(x,(,,fi 2 ) for Ji's 
distributions, which are integrated over the transverse momenta. 




Figure 1: Schematic diagram of the unintegrated generalised distribution H (x,£; k t , k' t ; /i 2 ) of 
([!]), showing the last step of the evolution. We use the 'symmetrized' variables of Ji 0. 



The cross sections are calculated from these universal H distributions using the k t fac- 
torization approach, and the explicit k t dependence allows the inclusion of the full kinematic 
behaviour. Assuming, for the moment, that k t has been appropriately integrated over, then, 
in the £ — > limit, the distributions reduce to the conventional diagonal quark and gluon 
distributions^ 



H q (x,0) 



q(x) for x > 
— q(— x) for x < 



(2) 



H g (x,0) = xg(x). 



There is a minor difference with the distributions of Ji in that H a = xHi 



The equations for the /x 2 evolution of the distributions of (0) are well known, but as the H a 
possess several arguments it is impossible to determine the starting distributions by fitting to 
the available sparse data. However here we show that it is possible to determine the x, £, k 2 , n 2 
dependence of the distributions from the values of the conventional integrated distributions, 
q(x,[i 2 ) and g(x,fi 2 ), that have been determined in the global analyses of deep inelastic and 
related hard scattering data, assuming that ^ 2 Cl and k' t ~ k t . For most practical applications 
indeed £ 2 <C 1, since it is specified by exp(— A77) where Ar] is the size of the rapidity gap. We 
also assume that \t\ <C /i 2 , where —t is the square of the momentum transfer. Indeed if \t\ were 
greater than the starting scale /i 2 ,, then the logarithmic structure of the evolution from /i 2 , to 
|t| would be destroyed. Although we cannot predict the t dependence, for small t it is expected 
that the t behaviour can be factored off as the proton form factor, see, for example, Ref . || . 

The procedure to determine generalised unintegrated distributions is to carefully combine 
two existing results. First, it has been shown that generalised (or skewed) distributions are 
completely determined from knowledge of the conventional distributions for small £ ||. Second, 
it has been found that the unintegrated distributions can be determined from the conventional 
(integrated) distributions -- the key observation is that angular ordering is only necessary in 
the last step of the evolution |5j . 

We briefly review these two results in Sections 2 and 3 below. In Section 4 we proceed to 
use this information to derive a formula for the unintegrated generalised distributions. Then 
in Section 5 we use the formula to compute the typical behaviour of the unintegrated skewed 
gluon distribution. To gain insight we also show other predictions for the gluon. Finally we 
give a simple phenomeno logical form that approximately reproduces the exact result. Section 6 
contains a brief conclusion. 

2 Generalised from conventional distributions 

Here we describe how the skewed (or generalised) distributions H a (x, £), with a = q or g, can be 
determined from the conventional parton distributions, q(x) or g(x), for £ 2 <C 1. For simplicity 
we omit the /i 2 argument. The essential result, due to Shuvaev 0, is 

H a (x, £) = £ dx' K a (x, £; x') f a (x') (3) 

with a = q or g, where the kernels K a are known and where, for £ <C 1, f a (x') reduces to the 
conventional distributions q(x') or g(x'). In this limit the Gegenbauer moments of the gener- 
alised distributions become equal to the Mellin moments of the usual parton distributions. For 
computational purposes, it is convenient to weaken the singularities in the integral expressions 
for K a by integration by parts. Then, for t = 0, (H) becomes 2 H 



H q (x,£) = / dx' — Im / 
J-i n Jo 



2 T r 1 ds 



* Jo y( s )^/l-y(s)x> 




(4) 



2 Note that here we adopt Ji's g] definition £ = (x\ — Xi), whereas in Ref. Q £ = (x\ — x%)/2. 

2 



H g {x,t) 



dx' 



where 



2 Im r l d8 ( x + t(i-*)) 

" Jo y(s)y/l-y(s)x' 
As(l-s) 



d f gjx'Y 

dx' \ \x'\ 



y{s) 



X 



-K(f-s)' 

For small x, where we may assume the conventional distributions behave as 

xf a (x) ~ N a x- X % 



(5) 



(6) 



(7) 



(|j) and (^) may be greatly simplified. Indeed we find that the ratios of the skewed to the 

diagonal distributions are 

H a (x,0 



where R a are simple known functions H. 



R a (x/£, A ) 



(8) 



3 Unintegrated from conventional distributions 



We outline the procedure || which enables the distributions f a (x,k 2 , \i ), unintegrated over 
k t , to be calculated from the conventional (integrated) distributions a(x,fi 2 ) 
xg(x,fi 2 ). We start from the DGLAP equation 



xq{x,fj, ) or 



da(x, /i 2 
<91n/i 2 



2n 



J2 [ A Paa'(z) a' (-,//) dz - a(x^ 2 ) V a (A) 



(9) 



where in the virtual contribution 



v%a) = y: 



1-A 



Pa'a(z) dz. 



(10) 



In the case of the g ^ gg splitting we have to insert a factor z in front of P gg in ([TJ]), to account 
for the identity of the produced gluons. The real emission term in (^j gives the unintegrated 
parton density, except that we must allow for the modification due to virtual effects. The 
virtual contribution does not change the parton k t and may be resummed to give the survival 
probability T a that parton a remains untouched in the evolution up to ft 2 . This survival factor 

is 

^ dk' 2 a s (k' 2 ) 

ik* 



T a (k tlf i) = expf-jT 



ki 



/2 



271 



V a (A) 



(11) 



and so the unintegrated distribution becomes 

da(x, k 2 ) 



/«(*,#, /O = T a {k t ^) 



T a (fc t ,/i) 



d In k? 

L J real 



(12) 



2vr 



E 



l-A 



Paa'(z) a \-,k t )dz 



X 



It is at the last step of the evolution, shown in ([12]) , that the unintegrated gluon becomes 
dependent on /i 2 . We have to make an appropriate choice of the cut-off A in ([TT|) and ([12]). We 
impose angular ordering in the last step of the evolution [|7|, § 



6(0-0') =► n > 



zkt 



\-z 
Thus the maximum allowed value of the integration variable z is 



i.e. A 



1 - z n 



H + kt 



(13) 



(14) 



/i + fct 

and similarly A = k' t /(fi + k' t ) in (PI). 

In summary, the /x dependence of the unintegrated distributions / a (x, fc 2 ,/i 2 ) enters at the 
last step, (|I~2"D, of the evolution through the survival factor T a , which results from the resum- 
mation of the virtual contributions, and through the cut-off A. 



4 Unintegrated generalised distributions 



To calculate the generalised (or skewed) unintegrated distributions H a (x, £; Af , /i ) from the 
conventional distributions a(x, /i 2 ), we combine the procedures outlined in Sections 2 and 3. 
First, we extend (f|) so as to describe the evolution of the skewed (integrated) distributions 
H a (x,^n 2 ) 



&ff a (x,£,/i 2 ) a s 



<91n/i 2 



2n 



i Jx \ X I \Z / Z 



(15) 

The splitting functions P aa i are given by eq. (22) of Ref. [Q], except that P gg (x,£) = xP gg (x,£), 
which arises because H g = xH^ 1 . For the antiquark contribution we should use —H q (—x/z, £, /i 2 ) 
in (^); see @. V\ and V2 are the virtual loop contributions corresponding, in the axial gauge, 
to the self-energies of the emitted and absorbed partons respectively, 



K = \ E r A% Pa>a(z)dz, 
2 ~, Jo 



(16) 



in analogy to ([10|). We will specify the cut-offs on the z integrations in a moment. 



Resumming the virtual contributions, as before, introduces the survival factors \fT[ and 
y/To~, where now a square root occurs due to the factor ~ in ([H]). In this way we obtain the k 2 
dependence of the skewed distributions, namely 



2 ,.2\ 



^a{x, q, k t , /! 



T a T a 



2vr 



E 



!- A dz 



& 



— P aa ,[ z ^—\ H a ,(-,£,ki 



with 



rpa 



cxp 



m 2 dkl a s (k' 2 ) 



2vr 



;i7) 



(18) 



Eq. (|17D shows the modification necessary to extend ([L2|) to skewed distributions. As in Sec- 
tion 3, the cut-offs in the z integrations of ( |17"D and (|16|) are chosen so as to respect angular 
ordering in the last step of the evolution. We obtain 



A 



A! 



1 + 



^n 



2x + f - £z n 



1 - z r , 



(19) 



(20) 



(1-^n 



1 + 



2£Zn 



2x + £ - 2£z n 



(21) 



where z max = fi/fi + k t of (P^|). Although we have assumed that the emitted and absorbed 
partons of Fig. 1 have the same kt, they have different momentum fractions, x\ ^ X2, and so 
the angular ordering gives different upper limits on the z integration in (17) , and in (|lfj) for V\ 
and\/ 2 §. 

We are now in a position to use the results of Section 2. We substitute (|]) and (|5p for the 
integrated skewed distributions H a i into the right-hand-side of flTT|). In this way we determine 
the unintegrated skewed distribution H a (x,£,; k^fi 2 ), for £ 2 <^ 1, in terms of the conventional 
distributions, q(x, /i 2 ) or g(x,fi 2 ). 



Note that there is no skewed effect coming from the z< 1 region of integration in (|17D . In 
this domain the second argument £z/x of the splitting function P aa r is very small, while the 
first argument (x/z ^> ^) of the integrated distribution H a i becomes large. As a consequence 
the £ dependence is negligible. Therefore we may replace the singular (1/z) part of P aa i by the 
(diagonal) BFKL kernel, and use this approach even for partons which incorporate the BFKL 
contribution in the small x region, just as was done in the diagonal case ||. Recall that in 
Ref. || it was seen that the unintegrated partons calculated from pure DGLAP partons (as 
in Section 3) are very similar to those calculated using partons which incorporate BFKL 
effects. It was concluded that the imposition of the angular-ordering constraint in the last step 
of the evolution is more important than the BFKL effects. 



5 Numerical results and discussion 



Typical behaviour of the unintegrated generalised gluon, H g , is shown by the continuous curves 
in Fig. 2. The results are obtained by evaluating (|P7|) for H g , using MRST99 conventional 
parton distributions |H| to first determine the integrated distributions, H a >, from @ and (H). 
We show results for the most relevant case, xi <C xi, that is for x = £/2. We take £ = 0.1, 0.01 
and 0.001 with fi = 10 GeV. To illustrate the /i dependence we also consider /i = 100 GeV for 
the choice £ = 0.01. In contrast to DGLAP k t ordering, the angular ordering, (|1^), imposed on 
the last step of evolution, enables the unintegrated distribution H g to extend into the k t > fi 



domain. However 
integral fll7D (or (|i 



H g decreases for large kt when z m3iX approaches x - - the lower limit of 
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Figure 2: The unintegrated skewed gluon distribution, H g (x,£ = 2x;/c^,/i 2 ), as a function of 
kf for various values of x and /x 2 . The choice £ = 2x corresponds to the physically relevant 
fractions x\ — £ and 22 = in Fig. 1. The continuous curves are the full result obtained from 
MRST99 partons |L0| by inserting (f|) and @ into (|T7|). The dashed, dot-dashed and dotted 
curves are shown only for information and correspond, respectively, to the values obtained from 
(1221), (E3I) and (M). 



To gain insight, we compare the results with different approximate formula for H g . The 
most naive possibility is to take 



H g (x,£;kt,ii 2 ) 



d(yg(y,k 2 )) 



d In k 2 



(22) 



y=x+£,/2 



This simplified distribution, which is widely used in the literature, is shown by the dot-dashed 
curves in Fig. 2. For not too small x and large k 2 , formula (^) gives a negative H g ; the reason 
is that the virtual term in Iffi) starts to dominate as x increases. This shortcoming may be 



avoided by developing a procedure originally implied in the work of DDT [|TTJ , see also [12 
For generalised partons, this leads to the formula 



Ha\ x i si k t , fi 



d 



d In kt 



T?T2a{y,k; 



In the diagonal case £ 



and Tf 



1 2 



Q T a 

dink? 



y=x+£/2 

T a . Then the derivative 



Z7l 



(23) 



(24) 



cancels the virtual contribution in (H), and H a is positive everywhere f|. However for the 
skewed case the situation is worse. Our sample curves correspond to x — £/2 for which A2 = 1 
and T 2 a = 1. Then only half of the negative term is cancelled by the derivative dy/Ti /dink 2 , 
and even the DDT- improved formula (p3|) may give negative H a ; see the dashed curves in Fig. 2 
for £ = 0.1 and k 2 > 100 GeV 2 . Note, however, the improvement of (^) as compared to ([22]) 
for the smaller £ values shown by Fig. 2. Due to the presence of the Tf survival factors in 
(|23|), the gluon distribution H g decreases for k t <C /i contrary to the behaviour of the simple 
derivative formula ([Z2"D. This reflects the small chance for low k t partons to remain untouched 
in the long evolution up to fi. 

In order to reproduce the behaviour of the unintegrated generalised gluon H g , for £ 2 <^ 1, 
in a simpler form than the full formula ([Ij]), h is useful to modify the DDT-like approximation. 
We retain the well-justified and important leading logarithm contribution coming from 



8VT9 
d In k 2 



N c o>s 
2vr 



In 



1 - z n 



(25) 



but tune the non-leading contributions and an additional valence quark contribution (needed 
for x ~ 0.1) to approximately reproduce the full result shown by the continuous curves in Fig. 2. 
The result is the phenomenological form 



H n 



„)£; k t , fj, 



R, 



dygjy, k 2 ) 2 



, N c a s A fi+jh 

+ yg{y,k f ) — — I ln — t j — + 1.2 



/// 



2tt 



h. 



^ + kl 



(26) 



+ 5 7T (y u ™\(y, k t) + yd va \{y, k 2 ) 



2tt 



y=i 



where R g is defined in (|j) and given in ||] . Form (|26|) gives the dotted curves in Fig. 2 and is a 
reasonable numerical approximation to the full result. It is a rough and ready way to estimate 
the values of H g for the physically relevant case x = £/2. 

Of course, from a formal point of view, we can add to the generalised distribution, H of 
(|17D , any function 

AH(x,0 = /(z,Oe(|-|*l) (27) 

with support only in the time-like ERBL region, |x| < £/2 ]T3| . Such a contribution is not 
obtained from the conventional partons determined from deep inelastic scattering (DIS) data. 
However for the processes mentioned in the Introduction we only need the unintegrated skewed 
distributions in the space-like (DIS) domain. 



6 Conclusions 

We have presented a procedure which allows the unintegrated skewed parton distributions, 
H a (x,£] fc 2 ,/x 2 ), to be calculated from the conventional distributions, g(x,/i 2 ) and g(x,n 2 ), for 
£ 2 <C 1 and any x. The key formula is fll7D , with (||) and (H) as input. A crucial advantage of 
using unintegrated distributions is that they enable the full kinematics of the hard subprocess 
to be retained, and in this way incorporate a major part of higher order QCD corrections. 
Moreover, the unintegrated skewed distributions are based on conventional partons fitted in 
global analyses of DIS and related data, and therefore have no uncertainties coming from the 
lack of knowledge of the (non-perturbative) initial conditions. They can therefore be used to 
make predictions for physical processes, including those for diffract ive production listed in the 
Introduction. 
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